The transition from adiabatic inspiral to geodesic plunge for a compact object around 

a massive Kerr black hole: Generic orbits 
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The inspiral of a stellar mass compact object falling into a massive Kerr black hole can be broken 
into three different regimes: An adiabatic inspiral phase, where the inspiral timescale is much larger 
than the orbital period; a late-time radial infall, which can be approximated as a plunging geodesic; 
and a regime where the body transitions from the inspiral to plunge. In earlier work, Ori and Thorne 
have outlined a method to compute the trajectory during this transition for a compact object in a 
circular, equatorial orbit. We generalize this technique to include inclination and eccentricity. 
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I. INTRODUCTION AND MOTIVATION 

Extreme mass ratio inspirals (EMRIs), in which stel- 
lar mass compact objects radiate gravitational energy 
and fall into their massive black hole companions, are 
promising sources of gravitational waves. LISA [1], the 
proposed space based gravitational wave detector should 
detect waves from the last stages of such inspirals. A 
clear theoretical understanding of the dynamics of EM- 
RIs is vital to the detection of these gravitational waves. 

The small mass ratios, which typically lie in the range 
/i/Af = 10"^ - 10~^, allow EMRIs to be treated within 
the framework of perturbation theory. The trajectory 
of the compact object can be roughly broken into three 
regimes: (a) An adiabatic inspiral phase, during which 
the dominant inspiral mechanism arises from the radia- 
tion reaction force on the smaller object. In this stage, 
the time scale over which the characteristic radial sep- 
aration (between the compact object and its black hole 
companion) changes is large compared to the orbital pe- 
riod. This allows us to approximate the trajectory as a 
sequence of bound geodesies, (b) A plunge phase, during 
which stable geodesies do not exist. It has been shown 
[2] that the effect of radiation reaction is negligible dur- 
ing the plunge and that this phase can be modeled as 
a geodesic infall. (c) A regime where the spiraling com- 
pact object transitions from adiabatic inspiral to geodesic 
plunge. The course of motion at this juncture shows as- 
pects of both, the self-force from radiation reaction and 
the effects of unstable geodesies. 

In [2], Ori and Thorne introduce a method to predict 
the motion when the object is constrained to an approxi- 
mately circular, equatorial orbit. We generalize this pro- 
cedure to include inclined and eccentric trajectories. A 
few modifications to the prescription in [2] are introduced 
to handle such generic orbits. The results from our gen- 
eralized prescription are in excellent agreement with [2]. 

The simple calculation described in this paper is meant 
to serve as a stopgap for many other open and important 
problems. There has been recent progress ([3, 4] and ref- 
erences therein) in the development of a code to solve the 



Teukolsky equation in the time-domain. The world line 
of the compact object serves as an input to this code. 
While the world line in the adiabatic phase can be calcu- 
lated from a frequency-domain based Teukolsky equation 
solver [5, 6], the trajectory in the transition regime for 
completely generic obits remains unknown. This calcu- 
lation will provide the missing link needed to generate a 
complete inspiral trajectory. 

A number of researchers are working towards solving 
the self- force problem exactly [7, 8]. Such an exact solu- 
tion can be separated (at least qualitatively) into time- 
reversal symmetric and asymmetric components. The 
symmetric component (the "conservative self-force" ) con- 
serves the integrals of motion. On the other hand, 
the asymmetric component (the "dissipative self- force") 
leads to non-zero time derivatives of the integrals of mo- 
tion. Recent advances demonstrate that we are mak- 
ing steady progress on this problem. For example, the 
self force is now essentially understood for circular orbits 
around Schwarzschild black holes [8] . Although approxi- 
mate, the results in this paper may serve as an indepen- 
dent check for these solutions. It is worth noting that if 
it becomes possible to include the conservative force in 
a simple way, we should be able to build its impact into 
the formalism developed here. This work may also be of 
interest for numerical relativity — a perturbative inspiral 
constructed by the techniques discussed here may be an 
accurate point of comparison for full numerical inspirals 
for small ratios (and may even be useful, if not so accu- 
rate, for mass ratios that are not strictly perturbative). 

Ref. [9] discusses the transition when the compact ob- 
ject is in an eccentric, equatorial orbit. However, the 
focus of that paper is to calculate the transit time and 
estimate the probability for LISA to observe such a tran- 
sition. Our intent is to generate the world line during 
the transition. We also choose our initial conditions dif- 
ferently than they are chosen in Ref. [9] ; we discuss these 
differences in more detail in Sec. III. 

The rest of the paper is organized as follows: Section II 
discusses circular orbits with arbitrary inclination. Sec. 
Ill generalizes the formalism developed in Sec. II to in- 
clude eccentricity. Finally, we summarize our results in 



2 



Sec. IV. 



II. THE TRANSITION TRAJECTORY FOR 
CIRCULAR ORBITS 



Up to initial conditions, a set of three constants, the en- 
ergy, E, the component of the angular momentum along 
the spin axis, L^, and the Carter constant, Q define a 
geodesic. The Carter constant has an approximate in- 
terpretation of being the square of the component of 
angular momentum perpendicular to the spin axis. As 
the compact object radiates, the "constants" that de- 
fine its geodesic will gradually evolve. (We will refer to 
[E{t), Lz(t),Q{t)] as the "constants", although they are 
slowly evolving.) A common approach to model the adi- 
abatic regime consists of treating the motion as the se- 
quence of geodesies [5, 6] defined by these evolving con- 
stants. As pointed out in [10], this limit amounts to a "ra- 
diative" or "dissipative" approximation. A true adiabatic 
approximation would be a sequence of orbits in which 
each orbit included conservative self corrections. Since 
we currently use purely geodesic orbits as our background 
motion (in lieu of a self-force enhanced description), we 
will refer to a sequence of geodesies as an "adiabatic in- 
spiral" throughout this paper. Thus, within the adiabatic 
approximation, the world line of a particle is computed 
by mapping [E{t), L,{t),Q{t)] to [r{t),e{t),(t>{t)]. The 
symbols r, 6 and (p are the usual Boyer-Lindquist coor- 
dinates. 

In contrast, the plunge can be treated as a single unsta- 
ble geodesic with almost constant E, and Q. Thus, 
the passage from adiabatic inspiral to geodesic plunge 
must contain both these features — slowly evolving "con- 
stants" and marginal stability. 



Kerr Geodesies 



The following system of first order equations describes 
geodesies in a Kerr [11, 12] geometry: 



dr 



(2.1) 
(2.2) 
(2.3) 

(2.4) 



The potentials can be expressed as: 
1 
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The parameters (r, 6, 0, t) are the Boyer-Lindquist coor- 
dinates, M is the black hole mass, n is the perturbing 
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2Mr + and a 



is the spin parameter of the black hole. The constants 

{E,Lz,Q) represent the actual energy, momentum and 
Carter constant (in units of M , M"^ and Af^ respectively), 
not the dimensionless versions of them. By introduc- 
ing the perturbing mass explicitly, our notation deviates 
from previous literature. We do this in order to show the 
dependence of the transition phase on the mass of the 
perturbing object. We also set G = c = 1 everywhere. 



B. The last stable orbit 

A standard but not unique definition of the "inclina- 
tion" of a Kerr geodesic is given by 



cosi- = 



Q = 



COS^ i 
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(2.9) 
(2.10) 



It is possible to use i to eliminate the Carter constant. 
Thus, any circular orbit can be parametrized by its radius 
(r) and inclination ((,). 

The last stable orbit (LSO) serves as an important ref- 
erence point — the inspiral is adiabatic well before the 
compact object crosses the LSO and is approximately a 
plunge well after the crossing. Since the transition occurs 
in the vicinity of the LSO, a preliminary step in our com- 
putation is to determine r and {E, L^, Q) at the LSO for 
a given inclination at the LSO, tLSO- Note that i changes 
with time because it is a function of [E{t), Lz{t), Q{t)]. 

Circular orbits satisfy 



i? = Oand 
dR 



R' = 



dr 
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(2.11) 

(2.12) 
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We must have R' = because the LSO Hes at an ex- 
tremuin of R. We also require that 

d'^R 



R" = 



>0, 



(2.13) 



for the cxtrcmum to be stable. This implies that the 
orbit will be marginally stable if R" = 0. Thus, the 
three equations R = R' = R" = can be numerically 
solved for a given ilso to yield r, E, and Q at the 
LSO. 



ing points, ^max and ^min, where 

< ^min < ^max < TT . 



(2.17) 



The potential problems posed by the turning points can 
be eliminated by reparametrizing 6. Following Ref. [5], 
we use 



z = cos 



Z- cos X 



(2.18) 



C. The constants in the transition regime 

We need a model of the phase space trajectory, 
[E{t), Lz{t), Q{t)\ near the LSO in order to compute the 
world line of the compact object as it transitions from 
inspiral to plunge. To this end, we Taylor expand about 
the LSO to obtain 

E{t) ~ £;lso + {t- tLSo)i^LSO , (2.14) 

L,{t) ~ L^.LSO + (t - tLS0)4,LS0 , (2.15) 
Q{t) ^ QLS0 + (i-iLS0)(<3LS0 + ^<3) 

+5Q , (2.16) 

which are natural generalizations of equations (3.4) and 
(3.5) of Ref. [2]. The overdot denotes differentiation with 
respect to t. We will later see that our initial condition 
for t amounts to choosing ti^so, the instant at which the 
compact object crosses the LSO. This choice is consistent 
with the procedure in Ref. [2] — Eq. (3.14) of Ref. [2] 
implies a choice of tLso — 0. 

The constant terms in Eq. (2.16), 5Q and 5Q. are 
needed to guarantee that the trajectory remains circu- 
lar as we enter the transition. As the notation suggests, 
these constants are small compared to Qlso and Qlso- 
They are discussed in more detail when we discuss initial 
conditions for the transition in Sec. II F. 

The expressions (2.14), (2.15) and (2.16) do not include 
conservative effects of the self force. Pound and Poisson 
[10] have demonstrated that this omission will lead to 
observationally significant changes. Inclusion of these ef- 
fects would effectively alter the potentials, Eq. (2.5) - 
Eq. (2.8) leading to slight deviations of {E, L2-Q)lso 
and tlso (for a given ilso) from their geodesic values. 
The exact impact of these effects will not be known until 
we know what the corrections arc. We will later see that 
our results posses all the expected qualitative features 
despite this handicap. Moreover, the prescription in [2] 
and its generalization presented here can easily incorpo- 
rate these effects once they are known. 

The fluxes at the LSO remain a parameter in our code. 
We use the code developed in [5] to provide us the di- 
mensionless fluxes, {M/^ifE, {M/fi^)L^ and {l/n^)Q at 
the LSO. Equivalcntly, we can use the expressions in [13] 
(with zero eccentricity) for the dimensionless fluxes. 
D. Reparametrization of the ^-equation 

Numerical integration of the ^?-equation warrants some 
care. The issue arises because d0/dt vanishes at the turn- 



where 



/?(.-.,)(.-._) = /?.--. 



~ 9 ' 



(2.19) 



and /? = a^(/i^ — E"^)/ fj?. The ^-equation of motion now 
becomes 



dx ^I3{z+ - z) 



dt 7 + a^Ez{x)/tJi ' 



where 



_ E 



(r2+a2)2 



2MraL^ 

A/Lt 



(2.20) 



(2.21) 



Equation (2.20) can now be integrated without turning 
points because % varies from to tt to 27r as varies from 



6'min to 6lmax and back to 6lmin- 



E. The prescription 

In keeping with our main objective of obtaining the 
world line [r{t), 6{t), (j){t)] through the transition regime, 
we eliminate r by dividing equation (2.1) by (2.4) and 
squaring the result to obtain 



R{r,x) 



F . 



Vt{r,xr 

One more time derivative gives the acceleration: 



(2.22) 



fr 
df^ 



d f R \ ^ d ( R \ dx/dt 



dr \Vj^J dx dr/dt 



(2.23) 



Ideally, Eq. (2.23) must have other additive terms pro- 
portional to non-zero powers of /i. This is analogous to 
Eq. (3.10) of [2]. Excluding this term amounts to ignor- 
ing the conservative self force. 

Since the transition phase is in the proximity of the 
LSO, we can Taylor expand F about tlsOj -E'lso, L^^lso 
and Qlso to obtain 
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F{r,L,,E,x,t) 



1 d^F 



6 dr^ 



(r - tlso)^ 



d^F 



LSO 



drdL, 



(1/2 - L2,LSo)(?' - r-LSo) 



LSO 



+ 



drdE 

Thus, the acceleration now becomes ^: 



^JE-E,so)ir-r,so)+-^ 



[Q - Q-LSo){r - tlso) • 



(2.24) 



LSO 



1 d'^F 



2 9r3 



(r - rLSo)^ 



d'^F 



LSO 



drdL, 



{Lz - iz,LS0) + 



LSO 



dx dr/dt 



drdE 



{E - -Elso) 



LSO 



(Q - Qlso) 



LSO 



(2.25) 



We have not expanded the second term in Eq. (2.23) 
because we do not know the value of x at r = Tlso a 
priori. Similarly, the ^equation takes the form 



# ^ V<l>{r, x) 
dt Vt{r,x) 



(2.26) 



The trajectory in the transition phase can now be com- 
puted by integrating equations (2.25), (2.26) and (2.20) 
from some starting point outside the LSO to some ending 
point inside the LSO, for a given <-lso, with time varying 
E, Lz and Q. 



where 



a = — 



To 



1 
2 

+ - 



1 a^ 

4af3 

a2 



R 



E2 



K{t) = - 



dLzdf 

5 d' 
dQdf 

1 dLz 
H/M df 

k|lso , 

(a/3Ko)"'/' 



LSO 

^1 + 



E d^ 



£ dEdr 



J LSO 

dLz/dt 
' {li/M){dT/dt) 



(2.30) 



(2.31) 

(2.32) 

(2.33) 
(2.34) 



Initial conditions 



The angles, (j) and x can be set to zero without loss 
of generality. Setting x = corresponds to starting the 
inspiral at ^ = ^min- 

The choice of initial radius depends explicitly on /x. In 
Ref. [2], the authors define parameters, a, /3, k, tq and 
Rq. These are used to scale out the perturbing mass from 
the equation of motion and initial conditions. Although 
we prefer to retain dimensions in the equations of motion, 
we specify initial conditions in a dimensionless form, in- 
dependent of ji. This will be useful be in interpreting our 
results and making comparisons with Ref. [2]. Following 
Ref. [2], we define 



X 



= (-) 

\m) 



2/5 



r - TLSO 

i?0 

-3/5 



(2.27) 

(2.28) 
(2.29) 

^ Note that Eq. (2.24) ignores terms of order (/i/M)^ and higher. 



T = 4^ 



\i/5 i-ii^so dr 



with f = r/M, i = t/M, E = E/fi, L^ = Lz/ifiM) and 
Q = Q/OiMf. 

These definitions reduce to those presented in Ref. [2] 
when t = 0. It is useful to observe that k does not scale 
with II. We evaluate dr/dt, a, /? and kq at 9 = 7r/2 — tLSO 
because we do not know 6'lso a priori. Notice that X and 
T are dimensionless. 

The smoothness of the transition implies that there is 
no fixed instant at which the transition starts or ends. 
Motivated by the choices in Ref. [2], we set T ~ —1 at 
t = and stop the numerical integrator when X < X^ = 
-5. 

In summary, our initial conditions are T = —1, = 
and X = at i = 2. Setting T = — 1 at t = allows us to 
calculate tLSO and hence E{0) and Lz(0) from equations 
(2.14) and (2.15). We then solve R{E,Lz,Q,r) = and 
dR/dr = to obtain r(0) and Q(0). This is analogous 
to Sec. IIIC of Ref. [2] where they enforce X = \/=T to 
determine X &i t = Q. 



It is important to keep |T| small enough that our Taylor expan- 
sion about the LSO remains a valid approximation. 
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The trajectory is adiabatic before the start of the tran- 
sition. At i = 0, we must impose the condition [14-16] 
that circular orbits remain circular even under adiabatic 
radiation reaction. Thus, requiring that R = dR/dt = 
and R' = d^R/drdt = leads to expressions (3.5) and 
(3.6) of [5] for r(0) and Q{0) respectively. 

We can now substitute (5(0) and Q{0) in Eq. (2.16) to 
obtain two independent equations, 

Q(0) = Qlso - tLsoiQhso + SQ) + SQ and{2.35) 
(3(0) - Qlso + SQ, (2.36) 

which can be used to evaluate SQ and SQ. 



G. Code algorithm and numerical results 

The previous sections developed the steps required to 
calculate the compact body's trajectory as it transitions 
from inspiral to plunge. We now summarize the algo- 
rithm that was actually used to implement this prescrip- 
tion: 

(1) Take iLSO as input. 

(2) Compute E and at the LSO. 

(3) Obtain E and Lz at the LSO from the code developed 
in [5]. We may also use the expressions in [13] (which re- 
duce to the results in [5] for circular orbits), which will 
be particularly useful when we generalize to eccentric or- 
bits. 

(4) Choose initial conditions T ~ — 1, = and x = 

(5) Calculate E{0) and ^^(0) from equations (2.14) and 
(2.15). 

(6) Solve for r(0) and Q{Q) by imposing R — and 
dR/dr = at t = 0. 

(7) Compute r(0) and Q(0) from equations (3.5) and 
(3.6) of [5]. 

(8) Substitute Q(0) and 0(0) in Eq. (2.16) to evaluate 
SQ and 6Q. 

(9) Use a Runge-Kutta integrator on (2.25), (2.26) and 
(2.20) to compute the coordinates at the next step. A 
time step of 6t ~ 0.05Af works well. 

(10) Update the "constants", Et+i = Et + E6t, L^.i+i — 
Lz,i + LzSt and Qi+i — Qi + {Q + SQ)5t. The subscript 
i denotes a discrete time instant. 

(11) Repeat steps (9)-(ll) until X{t) ~ -5. 

The primary objective of this calculation is to compute 
the world line of the compact object during the transi- 
tion. Figures 1 and 2 illustrate r, 9 and cf) motions of 
the compact object for a typical set of parameters. We 
also show a plunging geodesic matched to the end of the 
transition. 

Table I shows the parameters and transit times for a 
range of inclination angles. In general, we find that the 
transit time increases with inclination. However, the di- 
mensionless transit time AT remains approximately con- 



stant, 

AT -3.3 -3.4, (2.37) 

when Xe ~ —5 for all values of a and l. Again, this is 
a consistent generalization of the result in Ref. [2] where 
they find AT ~ 3.3 for all circular, equatorial orbits. 

Radial trajectory in the transition regime 
3.4| , , , , , , , , 1 



3.35- 



3.3 



3.25 



3.2- 



50 100 150 200 250 300 350 400 450 

t/M 

FIG. 1: Radial trajectory during the transition (black line) 
from inspiral to plunge for a compact object of mass jj, = 
10~^M in a nearly circular orbit around a black hole with 
spin a = 0.8M. The compact object crosses the LSO at time 
i-Lso ~ 137. 5M. The inclination of the orbit at tLSO is tLSO = 
37°. The red line is a plunging geodesic matched to the end 
of the transition. 



Angular trajectory In the transition regime 

2.5| 1 1 1 1 1 1 r 




50 100 150 200 250 300 350 400 450 




50 100 150 200 250 300 350 400 450 



t/M 

FIG. 2: Angular motion during the transition for a compact 
object around a spinning black hole with identical parameters 
as in Fig. 1. 
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TABLE I: Fluxes and transit times for different inclinations. We set a = 0.5M, /x = 10 ®M , M = 1, T(0) = -1 and Xe = -5. 





rLSo/M 




(M//x2)i,,Lso 


(l//x')(3z,LSO 


a 


/3 


i?0 




To 


t/M 


AT 




4.23 


-0.00457 


-0.0422 


-0.000572 


0.00311 


0.0327 


0.0699 


0.603 


2.80 


944.9 


3.36 


10 


4.26 


-0.00446 


-0.0409 


-0.00684 


0.00304 


0.0327 


0.0677 


0.604 


2.81 


952.4 


3.36 


20 


4.32 


-0.00415 


-0.0375 


-0.0241 


0.00284 


0.0329 


0.0615 


0.610 


2.82 


974.9 


3.36 


30 


4.43 


-0.00368 


-0.0323 


-0.0481 


0.00254 


0.0333 


0.0523 


0.618 


2.84 


1012.6 


3.36 


40 


4.59 


-0.00314 


-0.0262 


-0.0733 


0.00219 


0.0342 


0.0416 


0.630 


2.86 


1065.9 


3.35 


50 


4.78 


-0.002594 


-0.0198 


-0.0946 


0.00184 


0.0363 


0.0309 


0.643 


2.88 


1134.4 


3.35 


60 


5.01 


-0.00208 


-0.0139 


-0.108 


0.00152 


0.0403 


0.0211 


0.657 


2.90 


1217.9 


3.35 



TABLE II: Variation of transit time with perturbing mass, 
At/M. We set a = 0.9M, tLSO = 0.001°, M = 1, = -1 and 
= -5. Note that tlso = 2.32M. 





t/M 


AT 


10"^ 


118.9 


3.449 


10-* 


185.6 


3.397 


lo-'^ 


292.2 


3.375 


10"*^ 


461.9 


3.367 


10-^ 


731.3 


3.363 


10-8 


1158.6 


3.362 



H. Compeirison with Ref. [2] 

The results in Ref. [2] provide an important sanity 

check for the case of circular, equatorial orbits. However, 
we have to account for the minor differences between the 
two approaches. Ref. [2] makes the approximations 



(2.38) 
(2.39) 



which lead to 





d4> 


and 


dt ~ 


dt 


ISCO 


dr 


dr 




'dt ^ 


'dt 


ISCO 




dL^/di 



K = — 



{lj,/M){dT/dt) ' 

dLz/di\isco 
' {ti/M){dT/dt)isco 



(2.40) 



which is a dimcnsionlcss constant. In our prescription, 
dr/dt varies with time. This time dependence has to be 
enforced because dr/dt is a function of 9, whose value at 
the LSO is not known a priori. The circular, equatorial 
case in Ref. [2] does not suffer from this pathology be- 
cause = 7r/2 at all times. Thus, wc treat slowly 
varying function of time. Table II shows the transit times 
for a nearly equatorial orbit (iLSO = 0.001) and a range 
of mass ratios. As the mass ratio becomes smaller, the 
variation in k decreases, and the dimensionless transit 
time converges to the limit where k is constant. 

Our initial conditions differ slightly from those used 
in Ref. [2]. Effectively, they use the Taylor expansion of 



R{r) to solve dR/dr = and d^R/{drdt) = for r(0) 
and r(0) respectively. In contrast, we solve the equations 
exactly. This leads to differences of less than 1%. 



III. ECCENTRIC ORBITS 

The methods developed thus far only discTisscd circular 
orbits. We now extend this technique to include non-zero 
eccentricity. In the absence of radiation reaction, the 
geodesic equations admit bound eccentric orbits. These 
orbits are conventionally parametrized by the semi-latus 
rectum, p, and the eccentricity, e. The radial coordinate 
can now be expressed as 



r{t) 



P 



1 + ecos'ip{t) 



(3.1) 



The angle ip{t) is analogous to the eccentric anomaly and 
can be solved for numerically. The geodesic has turning 
points at -0 = O.tt. Deep in the adiabatic inspiral, the 
compact object's trajectory is well approximated by a 
sequence of orbits with slowly varying p{t) and e(t). 

Geodesies beyond the LSO do not have turning points 
(where dr/dt = 0). This changes the situation consider- 
ably because the parameters, p and e are not well-defined 
anymore. Thus, the trajectory ceases to have turning 
points somewhere during the transition from inspiral to 
plunge. We will later show that this feature is naturally 
buried in our model of the transition. 



A. The last stable orbit 

As with circular orbits, the last stable bound geodesic 
is an important reference in our procedure. The inner 
and outer turning points (rmin and rmax) of the LSO are 
related to clso and plso through 



r-min = — and 

1 + CLSO 

_ Plso 
1 - glso 



(3.2) 

(3.3) 



Our goal is to determine plso and the constants 
{E,LziQ) at the LSO for a given tLSO and clso- This 
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can be achieved by requiring that 

— = at r = rmin , (3.4) 
ar 

R = at r = Tmin and r = r^ax • (3.5) 

Recall that the function R is given by Eq. (2.5) and l is 
defined by Eq. (2.9). We require Eq. (3.4) to be satisfied 
because the inner most turning point corresponds to a 
local maximum of {—R)- Equation (3.5) enforces the 
compact object's velocity to vanish at the turning points. 
Equations (3.4), (3.5) and (2.9) can be solved numerically 
for p and {E,Lz,Q) at the LSO. Appendix A describes 
the details of this numerical procedure. 

B. The constants during the transition 

As with the circular case, our initial conditions are such 

that we effectively choose the LSO crossing to occur at 
t = ^LSO- This allows us to expand the constants about 
the LSO to obtain 

E{t) ~ £;lso + {t- iLSo)£^LSO , (3.6) 
Lz{t) — Lz,Lso + {t - t-Lso)Lz,i,so , (3.7) 
Q{t) ^ QLSO + (i-iLSo)QLSO- (3.8) 

Notice that we no longer need the corrections, SQ and 6Q 
because there are no additional symmetries to constrain 



Note that we only expand about the constants, not the 
r-coordinatc, because there is no unique r at the LSO. 
In the absence of the first three terms in Eq. (3.10), 
the equation of motion is simply a geodesic at the LSO. 
This is consistent with our intuitive notion of "expanding 
about the LSO" . The existence of turning points presents 
a complication while integrating Eq. (3.10) numerically. 
We present a method to tackle this in Appendix B. 

D. Initial conditions 

We need initial conditions for r and dr/dt before we 
start the numerical integrator. Motivated by the initial 
conditions for circular orbits, we set T ~ — 1 at f = 0. 



Q{0) and Q(0) ; E{t), Lz{t) and Q{t) are independent. 

As discussed in Sec. II C, equations (3.6), (3.7) and 
(3.8) do not include conservative effects of the self force. 
Just as the circular case, this will lead to a slight shift of 
(E, Lz, Q)hSO and plso (for a given clso and tLso) with 
respect to their geodesic values. Again, our motivation to 
stick with this approximation stems from the facts that: 
(a) These effects can be incorporated into our prescrip- 
tion once they are known, and (b) Our results show the 
generally expected behavior, at least qualitatively. 

Numerical methods to calculate the change in the 
Carter constant due to gravitational-wave backreaction 
have recently become available [17, 18]. Work is in 
progress implementing that result in the code we use to 
compute the rate of change of orbital constants [19]. For 
now, we use the approximate expressions for Q described 
in [13]; it will be a simple matter to update our code 
when more accurate Q results are available. 



C. The prescription for eccentric orbits 

Our next task is to derive equations of motion to map 
the phase space trajectory to an actual world line. The 
angular equations, Eq. (2.20) and Eq. (2.26), remain un- 
affected. Our strategy for the radial equation is to ex- 
pand the geodesic equation about (-Elso, -Z^z,lsOj Qlso)- 
This leaves us with 



(3.9) 



(3.10) 

I 

This amounts to choosing tLSO- We can now determine 
[E{0), Lz{0),Q{0)], which can be mapped to (p, e, i) at 
t = 0. This mapping is allowed because the trajectory is 
adiabatic before t = 0. The coordinates at any point on 
the geodesic defined by [E{0), Lz{0),Q{0)] can serve as 
our initial conditions. For simplicity, we choose 

r-(O) = , (3.11) 

^(0) = 0, (3.12) 

m = 0, (3.13) 
X(0) = 0. (3.14) 



fr 
df^ 
dF 
dr 



dF 

dr 



+ 



drdE 
dF 



dF dx/dt 
dx dr/dt _ 

{E - E'lso) + 



d^F 



LSO 



dr 



drdL 

{r, x; ■E'lso, -^z.LSO, Olso) 



d^F 



{Q - Qlso) 



LSO 
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The equations of motion can now be easily integrated 
across the LSO. 



E. Code implementation and numerical results 

Taking eccentricity into account changes our algorithm 
shghtly. We summarize the code's algorithm as follows: 

(1) Take iLSO f^^d blso a-s input. 

(2) Compute E, and Q at the LSO. 

(3) Obtain and Q at the LSO from the expressions 
in Rcf. [13]. 

(4) Choose initial conditions T ~ — 1, = and x = 
at i = 0. 

(5) Calculate E{Q), £^(0) and Q(0) from equations (3.6), 
(3.7) and (3.8). 

(6) Map [E{Q),LM,Qm to (p,e,i). 

(7) Set r = p/(l + e) and dr/dt = at i = 0. 

(8) Use a Runge-Kutta integrator on (3.9), (2.26) and 
(2.20) to compute the coordinates at the next step. A 
time step of 8t ~ 0.05M works well. 

(10) Update the "constants", Ei+i ~ Ei + Ei^soSt, 
L^^i+i = L^^i + L^^i^soSt and Qi+i = Qi + Qhso^t. The 
subscript i refers to a discrete time instant. 

(11) Repeat steps (9)-(ll) until X ~ -5. 

Recall that the local minimum of the potential R is 
less than zero for bound orbits and is greater than zero 
for a plunging geodesic. The minimum is exactly zero at 
the LSO. These conditions can be used as sanity checks 
while performing the numerical integration. 

Figures 3, 4 and 5 show a typical trajectory during the 
transition from inspiral to plunge. The compact object 
starts at the minimum of the last bound geodesic be- 
fore the plunge. The radial coordinate increases until it 
reaches a maximum where R = dr/dt = 0. Subsequently, 
it turns around and heads toward the minimum. After 
executing a number of "whirls" near the minimum, the 
trajectory becomes unstable, and thus plunges into the 
central black hole. The whirls are evident from the angu- 
lar trajectory plotted in Fig. 5. We also show a plunging 
geodesic matched to the end of the transition. Notice 
that the plunge spends quite a bit of time at r ^ 2.8M 
— much more time than the transition trajectory. This 
is because the radiation emission built into the transition 
trajectory's construction pushes it off this marginally sta- 
ble orbit rather quickly. 

Table III shows the various parameters and transit 
times for a range of eccentricities. Note that the pa- 
rameters a, (3, i?o, Ko and tq (which are defined in Sec. 
II F) are evaluated at plso- In general, we find that the 
transit time is proportional to a. This is not surprising 
because a is the first term in the Taylor expansion of the 
potential, R. We also observe some degree of correlation 
between the transit time and tq, the parameter used to 
define the dimensionless time. 



Radial trajectory in the transition regime 

11 1 , , , — — — , , ^ 




% 50 100 150 200 250 300 350 400 

t/M 

FIG. 3: Radial trajectory during the transition (black line) 
from inspiral to plunge for a compact object of mass /i = 
in an eccentric orbit around a black hole with spin 
a — 0.8M. The compact object crosses the LSO at time 
thso ~ 196. 7M. The inclination and eccentricity of the orbit 
at fLSO are ilso = 45° and clso ~ 0.6 respectively. The 
red line is an unstable geodesic matched to the end of the 
transition. 



Radial trajectory in the transition regime 

3.2 I , , , r- 



3.1 - 




^IsO 300 320 340 360 380 

m 



FIG. 4: Same as Fig. 3, but zooming in on the final "whirls". 



F. Comparison with Ref. [9] 

As mentioned in the introduction, there are differences 
between our generalized prescription and the method de- 
veloped in Ref. [9] , which only models the transition when 
the compact object is in an eccentric, equatorial orbit. 
First, we set our initial conditions at the start of the 
LSO, whereas Ref. [9] sets the initial conditions at the 
end of the LSO. This educated choice allows ReL [9] to 
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TABLE III; Fluxes and transit times for different eccentricities. We set a = 0.8M, = 10 ^, tLSO = 45°, M = 1, Ts = —1 and 

Xe = -5. 
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Angular trajectory in the transition regime Radiai trajectory in the transition regime 




FIG. 5: Angular trajectory during the transition for the same 
set of parameters as in Fig. 3. 



derive an analytic form for the trajectory. Second, we 
differ in the choice of final conditions. ^ 

In attempting to make comparisons with Ref. [9], we 
found a number of typographical errors. Thus, we extract 
the essence of the calculation in Ref. [9] and present it 
in a form that (hopefully) makes the errors obvious. We 
start by expressing the radial geodesic equation as 



+ V{r) = 0, 



(3.15) 



FIG. 6: Comparison of our trajectory with approximate an- 
alytic results from Ref. [9]. The compact object is in an ec- 
centric, equatorial trajectory with parameters eLSO = 0.6 and 
^ — 10-^M. Its mass is ^ — IQ-^M and is around a black 
hole with spin a — 0.8M. The black line shows our trajectory; 
the blue line is obtained from Ref. [9] . The observed deviation 
is because the approximation in Ref. [9] is somewhat more re- 
strictive than ours. 



The orbit is unstable if the local maximum of V{r) is 
negative. Define 



where 



S2 



(3.16) 



-Max{y(r)} = -y(w) 



(3.17) 



Note that this impHes T^'(rinax) = and V'(r,nax) < 0. 
3 See Sec. IID5 and Ref. [20] of Ref. [9] for a description their We Taylor expand Eq. (3.15) about the maximum of V{r) 
choice of final conditions. corresponding to some {E, Lz, Q) just beyond the LSO to 
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get 



djSr) 
dT 



+ V{r 

max 

)+6rV'( 

• max ) 



d{5r 
~ltt 



T2" 



the transit time is correlated with a, the cocfScicnt of the 
first term in the Taylor expansion of the radial potential. 

The code developed in [3] and [4] solves the Teukolsky 
equation in the time-domain and thus computes gravita- 
tional waveforms for almost any given trajectory of the 
,(3.18) compact object. We intend to generate waveforms by 
feeding the world lines calculated using this prescription 
to the time-domain Teukolsky equation-solver. The re- 
sulting waveforms will be useful for LISA data analysis 
routines. It is also possible to use these waveforms to es- 
(3.19) timate recoil velocities from mergers of compact objects 
with black holes. 



where 



6r = r{t) - w , (3.20) 

' max ' max 

= 2/|y"(r„,ax)| . (3.22) 
The solution of Eq. (3.19) in the regime of interest is 

r{t) = r^^ - VIts sinh ( - — —] , (3.23) 

where tc is an integration constant. We can compare 
our numerical solution with Eq. (3.23) by letting the two 
trajectories intersect at some arbitrary instant. This free- 
dom is equivalent to choosing initial conditions. For ex- 
ample, Fig. 6 shows the two trajectories near Tmax for 
which tc is choscin such that they intersect at t = 300M. 
The compact object has mass /U = 10~^M and is in 
an eccentric orbit with clso = 0.6 around a black hole 
with spin a = 0.8M. Notice that Eq. (3.23) is valid 
only in the immediate vicinity of fmax because {dt/dr) 
and {E,Lz,Q) are assumed constant. Inclusion of the 
time-dependence of (dt/dr) is crucial because it leads to 
time varying 7, which alters the natural timescale in Eq. 
(3.23). This explains the observed deviation at large val- 
ues of \Sr\. 



IV. SUMMARY AND FUTURE WORK 

The primary focus of this paper is to provide an ap- 
proximate model for the trajectory of a compact object 
as it transitions from an adiabatic inspiral to a geodesic 
plunge. We have presented a generalization of the pro- 
cedure in Ref. [2], where circular, equatorial orbits are 
treated. We derive approximate equations of motion [Eq. 
(2.25) and Eq. (3.10)] by Taylor expanding the geodesic 
equations about the LSO and subjecting them to evolv- 
ing E, Lz and Q. We can now readily integrate these 
equations numerically. Figures 1 and 2 show the radial 
and angular trajectories for a typical inclined, circular 
orbit. We also plot the plunging geodesic that it tran- 
sitions to. Figures 3 and 5 are analogous plots for an 
eccentric orbit. Our numerical experiments suggest that 
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APPENDIX A: THE LSO FOR ECCENTRIC 
ORBITS 

The following set of equations need to be solved in 
order to compute p and {E, L^, Q) at the LSO for a given 



inclination (t) and eccentricity, (e): 

R{r,E,L,) = 0, (Al) 

R(^^;^^,E,L}j = Oand (A2) 

^{r,E,L,) = 0. (A3) 



Recall that R is given by Eq. (2.5). The carter constant, 
Q can be eliminated using Eq. (2.9). Applying an it- 
erative technique to solve the above equations directly 
can lead to problems because the terms that do not con- 
tain r are identical in equations (Al) and (A2). We can 
skirt around this problem by solving the equivalent set 



of equations, 

Ri{r,E,L,) = R{r.E.L,)=0, (A4) 
R2{r,E,L,) = R(^-}-±^^^E,L)j-R^{r,E,L,) 

= Oand (A5) 
i?3(r,i;,L,) = ^{r,E,L,)=0, (A6) 



using the standard Newton-Raphson method described 
in [20] . This iterative procedure takes an initial guess for 
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the solution as input. We use 
Lzfi = rofcosi 



6.1(1 -a/2), 

1 - 2qv^ 



Eq 



l~2v^ + qv 



3i;2 

3 



2qv3 



and 



(A7) 
(A8) 

(A9) 



where q = a/M, r = \fMjr and = (ro, iz,o, -E'o)^ is 
our initial guess for S = {r,Lz,E)'^. Let Si denote the 
sohition at any given iteration. The algorithm consists 
of incrementing Si as follows: 



where 



6Si = 



Ji = 



Bi = 



Si+i = Si + Xx SSi 



Ji ^Bi , 

/ dRi/dr dRi/dE dRi/dL, 
dR2/dr dR2/dE dR2/dL, 
^dRi/dr dRa/dE dR^/dL, 

{ — Rl,i,—R2,i,—R3,i)'^ , 



(AlO) 
(All) 



(A12) 
(A13) 
(A14) 



and A ~ 0.1. The subscript "i" denotes that the expres- 
sions are evaluated at (ri,Lz,i,Ei). We stop iterating 



when \Bi\ < x, where a; ~ 10 The method outlined 
here works well for a large fraction of parameter space. 



APPENDIX B: NUMERICAL INTEGRATION 
ACROSS TURNING POINTS 

As mentioned in Sec. IIIC, Eq. (3.10) passes through 
turning points. The numerical integrator can accumu- 
late error when dr/dt — > 0. This section describes our 
algorithm to resolve the issue. 

Let tp denote the instant at which dr/dt = 0. The ra- 
dial motion is highly symmetric about the turning point. 
Thus, we must have, 



dr 
Iti 



dr 



(Bl) 



where e is an infinitesimal duration of time. When the 
radial velocity becomes very small, we exploit this sym- 
metry and set 



dr 
di 



tr, + St 



dr 
di 



(B2) 



tr,—St 



which is the discretized version of Eq. (Bl). 
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